Interplay of intra- and interband absorption in a disordered graphene 
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The absorption of heavily doped graphene in the terahertz (THz) and mid- infrared (MIR) spectral 
regions is considered taking into account both the elastic scattering due to finite-range disorder and 
the variations of concentration due to long-range disorder. Interplay between intra- and interband 
transitions is analyzed for the high-frequency regime of response, near the Pauli blocking threshold. 
The gate voltage and temperature dependencies of the absorption efficiency are calculated. It is 
demonstrated that for typical parameters, the smearing of the interband absorption edge is deter- 
mined by a unscreened part of long-range disorder while the intraband absorption is determined 
by finite-range scattering. The latter yields the spectral dependencies which deviate from those 
following from the Drude formula. The obtained dependencies are in good agreement with recent 
experimental results. The comparison of the results of our calculations with the experimental data 
provides a possibility to extract the disorder characteristics. 

PACS numbers: 72.80.Vp, 73.50.Bk, 78.67.Wj 
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I. INTRODUCTION 

Both gapless (massles) energy band structure and pe- 
culiarities of scattering mechanisms determine the re- 
sponse of graphene on the in-plane polarized radiation 
in the THz and MIR spectral regions. Such a response is 
caused by both the direct intersubband transitions and 
the intraband transitions of free carriers accompanied by 
scattering processes (the Drude mechanism), see the re- 
views [l|-y| and references therein. Contributions of these 
two mechanisms in heavily-doped graphene are described 
by the inter- and intraband dynamic conductivities, (Xintor 
and (Tintra, respectively, given by 



ReCTintor ~ '^^^^^ ~ ^^-P)' 



CTintraCl^) 



1 + iuiTp ' 



(1) 



where ep and Tp stand for the Fermi energy and the 
momentum relaxation time, ap — {e^ / ATi)epTp /fi is the 
static conductivity, and ^(e) is a unity step-like fuction. 
Here we have neglected a weak imaginary contribution 
Imcintcr [3] and the smearing of the interband absorp- 
tion edge associated with the finiteness of the tempera- 
ture, T, and the collisional damping of the electron spec- 
trum. When such a damping can be neglected at T = 0, 
the expression for RcCTintcr describes an abrupt jump at 
photon energy hu) = Isp due to the Pauli blocking ef- 
fect. The interband and intraband contributions are of 
the same order of magnitude at ujtp ^ \J \epTp j nfi — 1. 
Therefore, the overlap of these contributions disappears 
if epTp/h. > 20 (see Fig. la). Thus, the description 
of interplay between both contributions is necessary in 
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FIG. 1: (a) Spectral dependencies of Re(T(aj) described by 
simplified Eqs. (1) for different doping levels determined by 
eFfp/h. (b) Intra- and interband transitions between lin- 
ear branches of gapless energy spectra, (c) Screening of non- 
uniform charge in ungated graphene with field distribution 
shown by dashed curves, (d) Gated graphene structure with 
unscreened charge and distribution of weak field shown by 
dotted curves. 



the THz (or MIR, depending on the doping level) spec- 
tral region for the typical samples with h/epTp > 0.2. 
As shown below, the quasi-classical description of the 
intraband absorption and, therefore, the expression for 
fintra(w) bascd ou the Drude formula is not valid if huj 
exceeds temperature of carriers [1, @ . This is because of 
the variation of the relaxation time as a function of the 
electron energy over energy intervals ~ fiuj, where the 
intraband transitions take place (see Fig. lb). 

The constant value for the interband conductivity 
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RecTintor associatcd with the frequency-independent ab- 
sorption was observed for undoped suspended graphene 
in the visible range of the spectrum. [3] . By analogy 
with electrical transport in field-effect transistors, MIR 
spectroscopy allows for the control of the Pauli blocking 
effect using the electrical gating Q. In sufficiently heav- 
ily doped graphene, the Drude mechanism of absorption 
was observed near the interband threshold, up to the 
MIR spectral region, [gt-fllj . Despite of a series of theo- 
retical studies (see Refs. 12-14 and references therein), 
which have treated the problem of intra- and interband 
absorption at different levels of approximation, there are 
two drawbacks in the interpretation of the above-listed 
experimental data. First, the quantum description of the 
intraband absorption was not connected with the scatter- 
ing parameters determined from the conductivity mea- 
surements and, second, an inhomogeneous smearing of 
the interband threshold due to long-range disorder was 
not considered. Thus, a re-examination of the intra- and 
interband absorption processes in a disorded graphene is 
timely now. 

In this paper, we present a general consideration which 
is based on the Kubo formula and takes into account both 
the scattering due to finite-range disorder and a long- 
range unscreened variation of concentration. The latter 
contribution appears to be essential because the screen- 
ing of the potential spatial variations associated with the 
disorder with the in-plane scales longer than (or compa- 
rable to) the distance to gate is suppressed. This mech- 
anism is schematically illustrated in Figs. Ic and Id. As 
a result, the inhomogeneous mechanism of smearing of 
the interband absorption at huj ^ 2eF should be taken 
into account, together with the thermal smearing of the 
Fermi distribution and broadening of the spectral den- 
sity due to the elastic scattering on the finite-range dis- 
order. The dynamic conductivity in the high-frequency 
range (if uJTp > 1) is expressed via the Green's func- 
tion, which is averaged over the finite-range disorder, 
and the carrier distribution averaged over the long-range 
non-screened disorder. The averaged characteristics of 
Green's function were obtained using the phenomenolog- 
ical description of the static conductivity similarly to Ref. 
15. Considering the interplay of the intra- and interband 
contributions, we analyze the spectral, temperature, and 
gate voltage dependencies of the absorption efficiency. 
The obtained results are compared with the recent ex- 
perimental data. (10 , 11] . The fitting procedure demon- 
strates that the disorder characterization, including the 
evaluation of the parameters of high-energy scattering 
and long-range inhomogeneties, can be performed using 
the THz and MIR spectroscopy. 

The paper is organized as follows. In the next sec- 
tion (Sec. II), we present the basic equations which de- 
scribe the inter- and intraband mechanisms of absorp- 
tion and consider unscreened long-range disorder in gated 
graphene. In Sec. Ill we analyze the spectra of relative 
absorption and their dependencies on the temperature 
disorder and as well as on the doping conditions. The 



last section includes the discussion of experimental data 
and the approximations used, and the conclusions. 



II. BASIC EQUATIONS 

We start with consideration of the dynamic conductiv- 
ity based on the Kubo formula which takes into account 
the elastic scattering and the unscreened long-range dis- 
order. Since the absorption is determined by the real part 
of conductivity, Kea^ , we neglect the luiauj contribution. 
41 



A. Dynamic conductivity 

Absorption of the THz (MIR ) radiation propagating 
along the normal to the graphene layer is described by 
the real part of dynamic conductivity, [6j given by 



Reo-^ 



27re^ 



J dE Ue - /bh 



■huj 



(2) 



\{a\vAa')\H{E - eo:)5{E - 



hu) 



Here is the normalization area, /(£«) is the equilib- 
rium distribution over the states \a) with energies 
and |(q;|{;|| |q;')P is the matrix element of the in-plane 
velocity operator v\\ between the states a and a' con- 
nected by the energy conservation law described by the 
(5-function. The a-states are defined by the eigenvalue 
problem {h + V^)5'x"' = £a^x'^ expressed via the single- 
particle Hamiltonian h, a random potential Vx., and the 
two-component wave function ^/i^^ with I =1,2. The sta- 
tistical averaging of Eq. (2) over the random potential is 
denoted as (. . .). In order to perform this averaging, we 
introduce the spectral density function 



^,(^x,rx') = 5]5(e-£„)*L"'*|,"j;. 

a 

Thus, Kea^j given by Eq. (1) is transformed into 
27r(eu)2 



(3) 



Recr^ = 



L2 



■JdE{fE-fE+nu)Jd^iJd^2 (4) 



X ( tr 



As (Xi , X2 )(7|| Ae+Huj (x2 , Xi )o- 



where tr . . . means trace over the isospin variable, v = 
10^ cm/s is the characteristic velocity, and = {& ■ e) 
is written through the 2x2 Pauli matrix &, and the 
polarization ort e. 

Further, the spectral density function is expressed 
through the exact retarded and advanced Green's 
functions, ^^(xi,X2) and ^^(xi,X2), as follows 

Ae {xi,X2) = i t?| (xi,X2) - (xi,X2) /27r. These 
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Green's functions are governed by the standard equation 



[E-iO- v(a- ■ pi) - Uxi - t^xi] 
x^?| (xi, X2) = (xi - X2) 



(5) 



which is written through the potential energy with sepa- 
rated contributions from finite- and long-range disorder, 
Wxi and labeled as fr- and Ir- respectively. The aver- 
aging (• • •) in Eq. (4) should also be separated as (• • 



and 



I fr- 



it is convenient to introduce the variables 



(xi -|- X2)/2 — X and xi — X2 = Ax in Eqs. (4) and (5). 
Using these variables and neglecting a weak contribution 
of Vi)x we obtain i)x-i-Ax/2 — ^^x- As a result, after the 
replacement E ^ E + v one can separate the averaging 
over fr- and Zr-disorder contributions and Eq. (4) takes 
the form; 

27r f 

ReCT^ = —i^^f J '^^ i{fE+v)lr - {fE+v+h^)lr) (6) 
Ae (Xi , X2)<7l| Ae+Huj (x2 , Xi )(7l| 



X / iiAx<^tr 



fr 



where we take into accoimt that (• • depends only on 

I Ax| and use j dyL = . 

We consider below the high-frequency spectral re- 
gion, ujTp ^ 1, when the correlation function {■ ■ ■) fr 
can be factorizcd through the averaged spectral func- 



tions AaxE = (^£;(xi,X2)) . The averaging over /r- 

\ / fr 

disorder gives the contribution: 



/ 



dAx(tr[...])^ 



'EO-wApE+hu 



^ll)- 



Here the spectral densities ApE should be written 
through the averaged Green's functions, G^^ = GpE and 
GpE, which are given by the matrix expression: 



^pE 



E-v{a-Tp{] 



■^pE 



pE 



yR 

^pE 



(8) 



where the self-energy-function is determined in the Born 
approximation through the functions 



VpE 



dpi 

(27r;i)2'^|P-P^I 



(9) 



E+iO~vpi 
(P-Pl) f 1 



Pii ( 

01 V 



E+iO+vpi 
1 



E+iO — vpi E+iO+vpi 



Using the Fourier transformation of (uxiitx2)/r = 
exp[— (xi — X2)^/2Z/r], one obtains in (9) the Gaussian 
correlation function Wap expressed through the strength 
of potential and correlation length, u and Ifr- Because 
^pE depends only on the matrix {& ■ p), the averaged 



spectral density of Eq. (7) takes the form 

(a-p) 



a{+) , A-) 

^r>R T ^r 



^pB ^ ^pE 



2p 



{4e'-4e>) (10) 



1(±) 



27r \SpE T VpEP £pE ± VpEP 

where we introduced the renormalized dispersion law and 
velocity, Eps = E - and Vps = v{l + Vpe/'"p)^ re- 
spectively. 

For the averaging over Zr-disorder in Eq. (6) we intro- 

oo 

duce the population factor {fs+vjir = / dE' (j)E-E' Ie' 

— oo 

with the kernel (j)E-E' determined by the characteris- 
tics of disorder, see Sect. IIB. Thus, based on the 
two assumptions employed (smoothness of ir-potential, 
Vv 0, and high-frequency approach, ojtf ^ 1) we 
obtain Bea^ written through the triple integral: 



Re(7(, 



wL2 



j dE j dE'<f>E-E' ifE' - fE'+Hu.) 



^ ^2 i^^pE + "^pE^ (^pB+ftw + ^iJ+nS) ' (11) 



where the contributions oc ^p^^^pB^ ^^'^ *^ ^pB^^pB^ 
correspond to the intra- and interband absorption pro- 
cesses, respectively. Similar expressions for the homoge- 
neous case, without Zr-contribution, were considered in 
Refs. 13. 



B. Non-screened long-range disorder 

Here we evaluate the square-averaged potential of Ir- 
disorder which is caused by a built-in random poten- 
tial Wx- We take into account the screening effect in 
the gated structure with the graphene sheet placed at 
2; = on the substrate of the width d and the static 
dielectric permittivity e, see Fig. Id. The electrostatic 
potential V^iz is defined by the Poisson equation with 
the induced charge at 2 = determined by the con- 
centration of electrons in the heavily-doped graphene: 
«x = (4/^4) Ep ^(£B - £xp) = [(£F - v^)/V^vn]\ The 
dispersion law e^x = vp + is written here through 
the screened potential tix = Wx + V^2=o- Performing the 
Fourier transformation over x-plane one obtains the sec- 
ond order equations for and which correspond to 
the substrate, > z > —d, and to the upper half-space, 
2; > 0: 



d^ 



V> 

v< 



0, 



2; > 
0> z> -d 



(12) 



The jump of electric field at z = is defined by the charge 
distribution over graphene (z — > 0) 



dV> 



dz 



z=0 



dz 



z=0 



4 



= 47re- 



I 

1,2 



so that the potential remains homogeneous at 



(13) 



0: 



^qz=0 ~ ''qz=0 — ''qz=0 

z — —d is defined by the in-plane homogeneous back- 



= V^,_n = Vnr.=n- The bouudaiy condition at 



gate voltage Vg as follows: 



f\z— — d 



The last 



boundary condition is the requirement that the potential 
vanishes at z — > oo. 
The solution of this electrostatic problem can be writ- 
ten through exp(±(7z) and which depends on V^z^o 
through the concentration n^. Solving the system of the 
boundary conditions, we find that the screening potential 
at z = should satisfy the requirement: 



K{y) = 



Vq.=o = --^^dKiqd), 

sinhy 
y (sinh y + e cosh y) ' 



(14) 



where the dimensionless function K describes the gate- 
induced quenching of Vqz=o due to the fixed potential 
at z = —d. According to Eq. (13), 4'q is the nonlinear 
function of the screened potential written through and 
the electrostatic problem given by Eqs. (12) and (13) is 
transformed into the integral equation for Vqz=o deter- 
mined by (14). Further, we consider the case of heavily- 
doped graphene, when ep ^ |ux| and the electron-hole 
puddles are absent, so that ~ Hi?(l — 2vx/£f) is writ- 
ten through the averaged concentration np. As a result, 
the linearized dependency between ^'q and the screened 
potential Vq takes place 



47re np (5q,o 



2n 



q 

ep 



(15) 



Therefore, Eqs. (14) and (15) give the linear relation 



between Va 



Wq and the Fourier component of 



the screened potential v„ 



The solution for Vq takes form: 



l + 8{e/hv)'^epdK{qd)' 



(16) 



where screening is suppressed (i.e. Vq — Wq) if d — )■ 
and the complete screening (|uq| <C \wq\) takes place at 
d — >■ oo. 

The averaging over Zr-disorder in the distribution 
{fE+v)ir is performed after the Fourier transform of 
fp with the use of the relation {exp [iv^^r / h]) = 
exp [—{vr/hy /2], where v = \/Jv^)^. is the averaged 
Zr-potential. Under the inverse Fourier transformation 
(i.e. the integration over time) of the averaged distribu- 
tion one obtains the result {fE+v)ir = / dE' (j)p_Ei fpi 
which is written through the Gaussian kernel 



exp [-{^E/vf/2\ /{^/2^v) 



(17) 



with the half- width v. Using Eq. (16) one can write v 
through the averaged built-in potential Wq. We consider 




FIG. 2; Dimensionless averaged potential v versus d/U,. for 
diflerent qf =12.5 (1), 25 (2), 50 (3), and 100 (4). 



here a simple case of the Gaussian distribution of Wy, 
described by the correlation function 



i^qu/q'/;^ 



- L2^q+q^o27r;f,u;2 exp [-{qkrf /2] , (18) 



which is written through the correlation length lir and 
the averaged built-in potential w = \/ {w'^)ir- Using Eqs. 
(18) and (19) one obtains 



V 

w 



dxxe ^^/^ 



1/2 



[l+gp{d/kr)K{xd/kr)Y 



(19) 



and this ratio depends on d/lir and on the coupling con- 
stant, gp — {8e'^/hv)lir/Xp cx epkr, written through the 
Fermi wavelength Xp — vfi/ep. If d/Ur ^ 1 one ob- 
tains v/w — >■ 1 and in the region d/lir > 1 one obtains 
v/w ~ {\ + e) / gp which is independent on d/Ur- In 
Fig. 2 we plot the function (19) for the case of Si02 sub- 
strate if the coupling constant gp ~10 - 100 corresponds 
to the concentration range 4 x 10^^ - 4 x 10"'^^ cm~^ and 
if lir ^50 nm. With increasing of np or lij. at gp > 100, 
the relation v/W ~ ■\/(l + e)/gi? appears to be valid for 
the interval of d/lir > 0.1. 



III. RELATIVE ABSORPTION 

In this section we consider the absorption of radiation 
propagated along the normal to graphene sheet placed on 
semi-infinite substrate with the high-frequency dielectric 
permittivity n. Here we analyze the relative absorption 
coefficient 



167rRe(Ti, 



/^(1 + y^) c 



(20) 



neglecting the second-order corrections with respect to 
the parameter 4TT\auj\ / c^/k. Note that the reflection 
and transmission coefhcients, and Ti^ (these values 
are connected by the energy conservation requrement: 
Rlu + + £,uj = 1), involve the first-order correction 
oc AttIcTi^I/c^/k and the Recr^j contributions should be in- 
volved for description of Ri^ and Ti^. 
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FIG. 3: Spectral dependencies of relative absorption in homo- 
geneous graphene (normalized to Winter = ^^-i-oo) at tempera- 
tures 77 K and 300 K (solid and dashed curves, respectively) 

for concentrations 5 x 10^^ cm~^ (1) and lO" cm"^ (2). Ar- (1) and 10^^ cm- 
rows correspond to the doubled Fermi energies and dotted homogeneous case, 
curves are on u)~^ asymptotics for Drude absorption. 



A. Spectral dependencies 

First, we consider given by Eq. (20) after substi- 
tution the dynamic conductivity determined by Eqs. (9) 
- (11). Performing the multiple integrations in Retr^, 
one obtains the spectral, gate voltage, and temperature 
dependencies of relative absorption. For calculations of 
the renormalized energy spectra and velocity defined by 
Eq. (9) wc use here the typical parameters of /r-disorder 
corresponding to the sample with the maximal sheet re- 
sistance ^--^3.5 kr2 and the correlation length If^ ~7.5 nm, 
see 



FIG. 4: The same as in Fig. 3 at T =77 K taking into 
account inhomogeneous smearing of threshold due to non- 
screened variations of concentrations around 5 x 10^^ cm"^ 
^ (2). Dashed curves are plotted for the 
v=Q. 




15[ for details. Level of long-range disorder u, corre- 
spondent to lir ^ \f, is determined by the ratio d/kr, 
see Fig. 2; here we consider heavily-doped graphene be- 
cause a more complicated analysis is necessary for the 
case when electron-hole puddles are formed. [IQl 

The spectral dependencies of normalized to the 
high-frequency interband contribution, Winter, are plotted 
in Fig. 3 for the homogeneous graphene case, U = 0, at 
concentrations 5 x 10^^ cm"^ and 10^^ cm"^ for temper- 
atures T —77 K and 300 K. Similarly to the schematic 
spectra in Fig. la, one can separate contributions from 
intra- and interband transitions at huj ~100 meV and 
^150 meV for lower- and higher-doped samples with the 
minimal value of ^(^/^inter ^^-^0.15 and ~0.Q8, respectively. 
A visible deviation from the Drude spectral dependence 
(cx uj~'^) takes place due to quantum character of in- 
traband transitions if huj > T. Different signs of these 
deviations, which depend on huj, T, and np, appear due 
to complicated energy and momentum dependencies in 
Eqs. (9) and (10). The smearing of interband absorption 
around the Fermi energies (2e p are marked by arrows) is 
defined by both the scattering-induced broadening of the 
spectral density (10) and thermal effect, c.f. solid and 
dashed curves in Fig. 3. 

In addition, the inhomogeneous smearing of inter- 
band absorption threshold due to non-screened long- 
range variations of concentration is essential, as it is 
shown in Fig. 4 for T =77 K. Here we considered a strong 



100 200 300 
V p (meV) 



FIG. 5: Fitting of experimental data from Refs. 10 and 11 
marked as 1 and 2: (a) momentum relaxation rate Vp ver- 
sus energy vp and (b) sheet conductivity versus gate voltage. 
Squares and triangles are experimental points from Refs. 10 
(1) and 11 (2), respectively. 



Zr-disorder case with v ~30 meV for ni? = 5 x 10^"'^ cm"^ 
and V ~20 meV for np = 10^^ cm~^ (we used the de- 
pendencies of Fig. 2 with w ~50 meV and d/hr <0.1). 
Similarly to the thermally-induced effect, the unscreened 
disorder effect results both in smearing of the interband 
absorption threshold and in enhancement of the minimal 
absorption in the region between intra- and interband 
contributions. 



B. Comparison with experiment 

We turn now to comparison of our calculations with 
the recent experimental data, [13, Uli where both the 
spectral dependencies of relative absorption in large-area 
samples and the static conductivity a versus gate voltage 
Vg were measured. Using the phenomenological momen- 
tum relaxation rate Vp suggested in Ref. 15 and the ex- 
perimental data for hole conductivity from Refs. 10 and 
11, one can fit the dependency cr{Vg) as it is shown in 
Fig. 5. This approach gives us the /r-disorder scattering 
parameters for energies vp <200 meV and we can ex- 
trapolate these data up to ^400 meV energies, which are 
necessary for description of mid-IR response measured in 
Refs. 10 and 11. 

In Fig. 6 we plot the spectral dependencies of rela- 
tive absorption measured in Ref. 10 for graphene with 
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FIG. 6: Fitting of relative absorption spectra (normalized to 
Cintcr) to conditions of Ref. 10 at gate voltages -30 V (1), -60 V 
(2), and -80 V (3). Solid and dotted curves correspond to ho- 
mogeneous and inhomogeneous (with w ~150 meV) samples. 
Squares, triangles, and diamonds are experimental points for 
the cases 1, 2, and 3, respectively. Arrows correspond to the 
doubled Fermi energies. 



hole concentrations ^2.2, ~4.3, and ~5.8xl0^^ cm~^ at 
T =100 K. Solid curves are plotted for the case of ho- 
mogeneous sample with the scattering parameters deter- 
mined from the transport data shown in Fig. 5. These 
dependencies are in agreement with experiment both in 
the intraband absorption region (at Yiuj <150 meV) and 
above threshold of interband absorption, at ?iw > Isp. 
But the absorption in the intermediate region fuxi ~200 - 
300 meV appears to be suppressed in comparison to the 
experimental data. The minimal absorption increases if 
we take into account the unscreened Zr-disorder contribu- 
tion with V '-^SO meV, 75 meV, and 60 meV for curves 1, 
2, and 3, respectively. Such values oiv are realized if w is 
comparable to £p and djli^. < 0.1, i.e. micrometer-scale 
of inhomogeneities takes place, see Fig. 2. Thus, we have 
obtained a reasonable agreement with the experimental 
data. But the smearing of the interband threshold is 
stronger in comparison to experimental data because the 
Gaussian model does not describe a real Zr-distribution. 
A more accurate description is possible with the use of 
Ir-disorder parameters taken from additional structure 
measurements . 

Similar fitting of ^tj/Cintor measured in Ref. 11 for 
hole concentrations ^4.3, ~6.5, and ~8.6, and 11x10^^ 
cm~^ at T =300 K is presented in Fig. 7. We plot the 
spectral dependencies both for the homogeneous sample 
(solid curves) and for the inhomogeneous sample tak- 
ing into account the strong ^r-disorder (dotted curves) 
with V --120 meV, 105 meV, 90 meV, and 75 meV for 
curves 1, 2, 3, and 4, respectively. Now, good agree- 
ment with experimental data for interband absorption 
takes place (see the above discussion of Fig. 6 about the 
Zr-disorder contribution) while the intraband absorption 
appears to be 2 - 3 times stronger in comparison to the 
disorder-induced contributions calculated. An additional 
contribution may appear due to relaxation via emission 
of optical phonons of energy Hujq. Such a channel of in- 
traband absorption is allowed ii ep > and the rate 
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FIG. 7: The same as in Fig. 6 for conditions of Ref. 11 
at gate voltages -30 V (1), -60 V (2), -90 V (3), and -120 V 
(4). Squares, triangles, diamonds, and inverse triangles are 
experimental points for the cases 1, 2, 3, and 4, respectively. 



of transitions is proportional to the density of states at 
£f — hujQ. As a result, this process may be essential for 
the conditions of Ref. 11, where >250 meV, while this 
mechanism should be negligible at lower concentrations 
used in Ref. 10. More detailed experimental data, for 
the spectral range 200 - 400 meV at temperatures <300 
K, are necessary in order to verify this contribution. 

Overall, the comparison performed here demonstrates 
good agreement with experimental data, which opens a 
way for verification of scattering mechanisms and long- 
range disorder parameters. A complete verification can 
be performed under complex analysis of transport (con- 
ductivity versus Vg and T) and spectroscopic data to- 
gether with a structural (e. g., STM) measurements of 
long-range disorder. 



IV. CONCLUSIONS 

We have examined the spectral, concentration (gate 
voltage), and temperature dependencies of the relative 
absorption due to intra- and interband transitions in 
graphene taking into account both elastic scattering 
and non-screened long-range variations of concentration. 
Smearing of the interband absorption edge, deviation 
from the Drude spectral dependence due to quantum 
regime of intraband transitions, and interplay of these 
two contributions are demonstrated. The results are in 
agreement with the recent experimental data obtained 
for heavily-doped graphene. [ifl, [Hi This analysis allows 
to extract disorder parameters in the high-energy region 
(up to huj above the Fermi energy). 

Let us discuss the assumptions used in the calculations 
performed. The main restriction of the results is due to 
consideration of the high-frequency spectral region, when 
Keauj can be written through the matrix product of the 
averaged spectral density function (10). More compli- 
cate description, based on the Bethe-Salpeter equation, 
is necessary for the low-frequency region, which is be- 
yond of our consideration. We also restrict ourselves by 
the phenomenological models for the finite-range disorder 
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scattering and for the long-range variations of concentra- 
tion. It is enough for description of spectral, gate volt- 
age, and temperature dependencies of £,uj because they 
are expressed through simple correlation functions which 
are similar for any microscopic nature of disorder. In ad- 
dition, a long-range random strain (l7l | may affect on the 
interband transitions in large-area samples. The scat- 
tering processes caused by phonons and carrier-carrier 
collisions should be analyzed under calculations of the 
intra- and interband absorption in clean samples. A pos- 
sible contribution of optical phonons is discussed in Sect. 
IIIB (more detailed analysis is beyond of the scope this 
paper) while the interaction with acoustic phonons is neg- 
ligible. The comparison with experiment was performed 
here for the hole doping case and electron-hole asymme- 
try reported in Refs. 10 and 11 is unclear. According 
to Refs. 18, the Coulomb-induced renormalization of re- 
sponse and scattering processes give a weak contribution 
for typical disorder levels. This disagreement requires an 
additional investigation. 



To conclude, we believe that our results open a way for 
characterization of scattering processes and long-range 
inhomogeneities using THz and MIR spectroscopy. More 
important, that effect of disorder on THz and MIR re- 
sponse is the main factor which determines characteris- 
tics of different devices (e.g., lasers [l^ and photodetec- 
tors [20|) in this spectral region. We believe that our 
study will stimulate a further investigation of these de- 
vice applications. 
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